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Abstract 

Arm locking is a proposed laser frequency stabilization technique for the Laser Interferometer Space 
Antenna (LISA), a gravitational- wave observatory sensitive in the milliHertz frequency band. Arm locking 
takes advantage of the geometric stability of the triangular constellation of three spacecraft that compose 
LISA to provide a frequency reference with a stability in the LISA measurement band that exceeds 
that available from a standard reference such as an optical cavity or molecular absorption line. We 
have implemented a time-domain simulation of a Kalman-fllter-based arm locking system that includes 
the expected limiting noise sources as well as the effects of imperfect a priori knowledge of constellation 
geometry on which the design is a based. We use the simulation to study aspects of the system performance 
that are difficult to capture in a steady-state frequency domain analysis such as frequency pulling of the 
master laser due to errors in estimates of heterodyne frequency. We find that our implementation meets 
requirements on both the noise and dynamic range of the laser frequency with acceptable tolerances and 
that the design is sufficiently insensitive to errors in the estimated constellation geometry that the required 
performance can be maintained for the longest continuous measurement intervals expected for the LISA 
mission. 

PACS numbers: 04.80.Nn, 95.55.Ym, 07.05.Dz, 07.05.Tp 
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I. INTRODUCTION 

The Laser Interferometer Space Antenna jl, [2| is a planned facility for observing gravitational 
radiation in the milliHertz frequency band, a regime rich in astrophysical sources. The LISA 
measurement concept 3] calls for laser interferometry to be used to measure fluctuations in the 
distance between freely-falling test masses contained within spacecraft separated by ~ 5 x 10 9 m 
with a precision of ~ 10 x 10 _12 m, or ~ 10 pm. The interferometric measurements are performed 
as a series of one-way measurements between pairs of spacecraft (SC) and then combined using a 



,1 



to form observables that suppress laser 



technique known as Time Delay Interferometry (TDI) 3, 
frequency noise while retaining gravitational wave signals. 

The capability of TDI to reject laser frequency noise is chiefly limited by imperfect knowledge 
of the absolute light travel times between the spacecraft (often referred to as the "arm lengths"), 
which is expected to have an accuracy of ~ 1 m or ~ 3 ns. With this level of arm length accuracy, 
the contribution of laser frequency noise in the TDI observables will satisfy the allocated equivalent 
path length noise of 2.5 pm/ a/Hz so long as the input laser frequency noise does not exceed a level 
of 



Hz \ / /2.8mHz\ 4 

Vpre-TDl(f) = ( 282 — ) ■ \j 1 + ( ) , (1) 



/ w V V / 

where the Fourier frequency / ranges over the LISA measurement band, 0.1 mHz < / < 0.1 Hz. 
The expected free-running noise level of the LISA lasers in the measurement band is roughly 
lOkHz/v^Hz ■ (1Hz//), which exceeds the requirement in (CQ) by more than four orders of mag- 
nitude at the low end of the LISA band. Consequently, the lasers must be stabilized using an 
external frequency reference. A number of candidate stabilization schemes have been studied and 
determined to be viable from a noise performance perspective (5, 6|. The current focus is on eval- 
uating other aspects of each candidate scheme such as complexity, robustness to implementation 
errors, and operational constraints so that the most effective design can be selected. 
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Two of the candidate schemes rely on arm locking, which utilizes the existing LISA science 
signals to derive a frequency reference from the geometry of the constellation. In one scheme, 
arm locking is the sole method employed to stabilize the laser frequency where in the other it is 
combined with another stabilization method in a hybrid system. We focus on the latter case in 
this paper. 

Since its original introdnctionfl, the a.™ .oc kl „ g concept has been rennedfl Q .eading to 
improvements in its exp ected performance. Arm locking has also been studied in a number of 
hardware models 
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12|, which have helped to identify potential implementation issues that were 



not readily apparent from frequency- domain studies of arm locking. Chief among these was the 
discovery that the inability to predict the heterodyne frequencies of the LISA science signals due 
to imperfect knowledge of the inter-spacecraft Doppler shifts leads to "frequency pulling" of the 
arm locked laser. If not properly mitigated, this frequency pulling can be so severe that the laser 
exceeds its dynamic range in a matter of hours. 

In this paper we present the results from a series of time-domain simulations of arm locking. The 
goal is to combine the attractive features of the frequency domain models (realistic noise sources, 
orbit models, etc.) with those of the hardware models (sensitivity to transients, non-linearities, 
etc.). The rest of the paper is organized as follows. In section [III we describe the problem of arm 
locking in LISA, defining the relevant signals. In section 1X111 we briefly review the frequency-pulling 
effect. We discuss our particular arm locking design and the technical details of our simulation in 
section |IV] and present results in section [V] 



II. ARM LOCKING MODEL 



To maintain the focus on the arm locking dynamics, we make a few simplifications in our model 
of the LISA interferometry. The first is that we only consider the interferometric measurements 
made between different SC (the "long-arm" signals) and ignore the additional measurements made 
between the SC and the proof mass (the "short-arm" signals). While a combination of both signals 
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are needed to reach the ~ 10 pm sensitivity levels required for detecting gravitational waves, the 
~ nm sensitivity level of the long-arm signals is more than sufficient for frequency stabilization at 
the level of (□). 

A second simplification is that we model the two lasers on board the master SC as a single laser. 
In reality one of the lasers will be phase locked to the other using a high-gain phase lock loop. The 
residual noise in this phase lock loop is expected to be far below the other noise sources considered 
in this paper. Readers interested in additional detail on the LISA interferometric measurement 
concept should consult one of the many overview papers 



A. Notation 

Our notation system is an adaptation of that used by McKenzie, et al.Q. The most notable 
difference is that we represent signals as fluctuations in frequency rather than fluctuations in phase 
and replace with v to reflect this. The three LISA spacecraft (SC) are labeled SCi, i = 1,2,3 
and it is assumed that SC\ is the master SC. The many different frequency signals are labeled 
with both an alphabetic and a numeric subscript. The alphabetic subscript refers to the physical 
nature of the signal while the numeric subscript refers to the SC involved in producing the signal. 
In a two digit numeric subscript, the first digit indicates the receiving SC while the second refers 
to the transmitting SC. For example, vsiz denotes the shot noise on the photoreceiver on board 
SC\ that is receiving signals from SC3. Appendix [A] contains a table summarizing the notation 
used in this paper and, where possible, the corresponding notation in (9)]. 



B. A single LISA laser link 

Figure [1] shows a schematic of the LISA constellation. We begin our analysis of the arm locking 
signal chain with Laser 1 on SC\, which produces light with a frequency vq\- This frequency is 
a combination of the intrinsic frequency noise of the laser, u^i, and the control signal provided 
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by the arm locking loop. As the laser departs SC± in the direction of SC3, it picks up a Doppler 
shift due to the motion of SC\. The magnitude of this Doppler shift is A" 1 V^i ■ 1713, where A 
is the wavelength of the laser, V^i is the velocity of SCi, and f) 13 is the unit vector along the 
path from SC± to SC3. For the purposes of calculating Doppler shifts, we make the assumption 
that fjij = —fjji even though the rotation of the constellation causes these angles to differ on the 
~ /iradian level 
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FIG. 1: Schematic of frequency signals relevant to arm locking. See section|TT]of the text for details 
and Appendix lAl for a key to notation. Adapted from Figure 1 of 



The laser then experiences a delay of 713 011 the order of 5 x 10 9 m/c ~ 17 s as it travels to SC3. 
At SC3, the signal picks up another Doppler term due to the motion of SC3. At the photoreceiver 
on SC3 it is interfered with the local laser (with frequency U03) to generate an electrical heterodyne 
signal with frequency 
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(2) 



where ^531 is a shot noise contribution due to the low light level of the received beam. A device 
we will refer to as the "frequency meter" (although it is more commonly called a phase meter Jl^ ) 
is used to measure the frequency of the heterodyne signal. The first step is digitization, which 
introduces a noise term due to fluctuations in the frequency of the oscillator used to drive the 
digitizers. To first order, this clock noise is additive with a spectral density that is proportional to 
the instantaneous heterodyne frequency, 



Here 2/3 (/) represents the spectrum of fractional frequency fluctuations of the clock on board 
SC3. The frequency meter also measures the frequency of the heterodyne signal relative to some 
model signal ^M3i- The model signal can be used to remove the slow drift of the heterodyne 
frequency caused by time-varying Doppler shifts or to impose a constant frequency offset in a 
phase-lock loop. The output of the frequency meter is given by 



C. Doppler Shifts 

The LISA SC will experience relative velocities along their lines of sight of several meters per 
second, resulting in Doppler shifts of several MHz. These Doppler shifts are approximately constant 
in the LISA measurement band and it is convenient to remove them prior to implementing arm 
locking. We begin by separating the SC velocity terms into a deterministic term arising from the 
SC orbits (V Oi) and a stochastic term arising from attitude jitter of the SC (SVi) 



£c3i(/) = "mi ■ j/s(/). 



(3) 



fA3l(*) = VHSl{t) ~ fj\/3l(*) + fC3l(*)- 



(4) 
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f i (t) = foi(t) + 8f i (t). (5) 

This in turn leads to two Doppler contributions in the heterodyne signals Umj, an orbital motion 
term {voij) and a SC jitter term iyjij) given by 




(6) 



vjij(t) = A 



sfi(t) - 5\? j(t - t^)] -fjij- 



(?) 



With these signal definitions, the heterodyne signal on board SC 3 can be written 



vmi{t) = voz{t) - v \(t - Ti 3 + v mi {t) + V J31 (t) + U S 3i(t) 



(8) 



D. Phase locking on SC2 and SC3 

Arm locking requires that the slave SC (SC2 and SC3 in our example) operate in a transponder 
mode, returning to the master SC a copy of the light field that was received. This is accomplished 
by using a phase-lock loop to control the lasers on board the slave SC. Using SC 3 in Figure [1] 
as an example, the controller G 3 adjusts the frequency of Laser 3 to minimize the output of the 
frequency meter, 1^431- In the Laplace domain, the output of Laser 3 will be 



"03 (s) = 



G 3 



VD3l{s) - Vj3l(s) - ZAj3l(s) - ^C3l(s)] + 



1 



1 + G3 



1 + G3 



(9) 



Under the assumption of a high-bandwidth phase lock loop, G3 3> 1, this simplifies to 



u Q3 (s) « z/ M3 i + voi(s)e 



(10) 
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E. Formation of the arm locking error signal 



The signal from Laser 3 is transmitted back to SC\, picking up a Doppler contribution form 
SC 3 , a time delay r 31 , and a Doppler contribution from SCi. At SC\ it is interfered with Laser 1 
on a photoreceiver, generating shot noise i>si 3 . Fluctuations in the heterodyne signal are measured 
by a frequency meter, which subtracts a model VMi^{t) and adds a clock noise vq\ 3 to produce the 
main science signal for the SC\ — SC 3 arm, ^413 (t). Using ([TO]) to replace uq 3 , Vais can De written 
as 

VA13{t) = [l>Ol(t) ~ V0l(t ~ 7i3 - 73l)] + [l/ji 3 (t) + Vj3i(t ~ T 3 l)] 
+ [fS13(*) + (< - T-3l)] + [fC13(*) + fC3l(* - 7-3l)] 

+ [fU13 (<) + ^D3l(i -7-3l)] - [fJlfl3(*) +VM3l(t -T31)]. (11) 

The second to last bracketed term in (fTTj) represents the deterministic part of the heterodyne 
signal. The model signal in the frequency meter, ^mi 3 (£) 5 can be used to remove this term, leaving 
behind a (hopefully small) residual error term, VEuit), 

VMn(t) = V m3 (t) + U D31 (t - T 3 i) - U M 3l{t - 73l) - ^13^)- (12) 

These residual errors lead to frequency pulling of the master laser. Section [TTTJpresents estimates 
for the size of these errors. With the deterministic terms (mostly) removed, vaiz can De represented 
in the Laplace domain as 

^avs(s) = u i(s) [l - e -'(™+™)] + [ UJ13 ( 8 ) + vj 3 i(s)e- s ^} + [is ci3 (s) + i/ C31 (s)e-"™] 

+ h?i 3 (s) + v sn (s)e~ s ™] + u E13 (s). (13) 

The SC\ — SC2 arm produces a signal, Vai2(s), that is analogous to z/413 (s). The arm locking 
sensor is a linear combination of these two signals that is used to estimate the Laser 1 fluctuations, 
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Voi, so that they can be suppressed in a feedback loop. The output of the arm locking sensor, 
labeled Ugi in Figure [1] is given by 



vayi 

VA13 



(14) 



where S is the arm locking sensor vector that describes the specific linear combination of the two 



individual arm signals. For example, the "common-arm" sensor, uses the sensor vector S_ 



Table 
literature 



I in Q 



L2> 2) 



provides expressions for several arm locking sensors that have been studied in the 



F. Noise Levels 



1. Intrinsic Laser Frequency noise 



For this work we assume a pre-stabilized, frequency tunable laser source with a frequency noise 
spectral density in the LISA band of 



u L {f)= (800^|) ■Jl + ( 2 ' 8 p 4 0.1mHz</<0.1Hz. (15) 



This is representative of the noise-floor of the Mach-Zehnder interferometer stabilization 
15 1 that will fly on LISA Pathfinder [l(3|, a LISA technology demonstrator mission. 



system 



2. Shot Noise 

Shot noise is uncorrelated at each detector and has an equivalent frequency noise spectrum of 



(16) 



XPrec VVlfe/ V 1Hz 
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where A = 1064 nm is the laser wavelength and P rec ~ 100 pW is the received power. For these 
numbers, (fT6l) gives Vg = 43 /^Hz/v^Hz ■ (//l Hz). 

3. Clock Noise 

The spectral density of the fractional frequency variations of the SC clocks are estimated to be 
(Table IV of Q) 



y(f) = 2.4 x lO-^/V?- (17) 
While LISA will employ a clock-transfer scheme to correct for differential clock noise between 



the SC 



17], we assume that that correction takes place in post processing on the ground and is not 



applied to the arm locking error signals on board the SC. 
4- Spacecraft Jitter 

The spacecraft jitter noise model is based on simulations of the drag- free control performance 
of LISA[18|. The jitter can be divided into two orthogonal components in the plane of the LISA 
constellation that are independent. Each of these has a position jitter in the LISA measurement 
band of 

6x(f) = 2.5nm/Vhfe 0.1 mHz < / < 0.1 Hz. (18) 

Note that due to the fact that the interior angle of the constellation is not 90 deg, the spacecraft 
jitter contributions from SCi will be partially correlated in the frequency meter signals 1^412 (^) and 
v ais {t)- Finally, we note that it would in principle be possible to remove the spacecraft jitter by 
including the "short-arm" interferometers in both the phase-lock error signals in the far SC and the 
arm locking error signals in the master SC. This would reduce the jitter from ~ nm to ~ pm in the 
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LISA band. As with the clock noise correction, this would require additional on-board processing 
and is not necessary to reach the pre-TDI noise requirement specified in <^j. 



III. LASER FREQUENCY PULLING AND HETERODYNE ESTIMATION 

In section [Til we explained how Doppler shifts arising from the SC orbits enter the long-arm 
frequency meter signals and how the deterministic parts of the signals are removed using models of 
the heterodyne frequency. If we take the expression for the main science signal of the SC\ — SC3 
arm, 1/413, as expressed in (fl3|) and consider only the terms resulting from the laser frequency, Uqi, 
and the errors in the heterodyne estimate, ueu, the result is 

^13(5) = 1/01(5) [1 - e-C*-™-™)] + ^3(5). 

If we take the low frequency limit, s — > 0, we find that the first term vanishes. In other words, 
the signal in uau is insensitive to fluctuations in laser frequency at zero frequency. The second 
term, however, is unaffected. The situation is obviously the same for vau and also for any arm 
locking sensor formed as a linear combination of VA12 and Vaiz- If the arm locking controller has 
any gain at zero frequency, it will cause the laser to ramp in an attempt to zero out the heterodyne 
error terms. This laser frequency pulling can be mitigated by reducing the arm locking loop gain 
below the LISA measurement band, although care must be taken to ensure that sufficient gain 
is still present within band. The rate of pulling is proportional to the error in the estimate of 
the heterodyne frequency, hence the design requirements of the control filter will be driven by the 
accuracy with which the heterodyne frequency can be estimated. 

As pointed out by ij, it is convenient to combine the heterodyne models from the individual 
arms into a common and differential component. For the case where SC\ is the master SC, the 
common and differential heterodyne models are 

u M +{t) = v Mi2 (t) + u M13 (t) 

m-(t) = VM12{t) - VMl${t) (19) 
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It is also expected that arm locking may require periodic re-acquisition, either because of an 
external disturbance (e.g. pointing of the high-gain antenna) or because some component of the 
arm locking system is in a non-desirable operating range (e.g. laser near a longitudinal mode 
transition, arms close to equal, etc.). Consequently, the heterodyne frequency only needs to be 
estimated for periods on the order of weeks. For such periods, it is appropriate to use a quadratic 
model, 

VMx{t) ~ "ox + loxt + a 0x t 2 , x = (+, -) (20) 

A. Expected Doppler 

Although the models for the heterodyne signals include both Doppler shifts and intentional 
frequency offsets, the Doppler shifts provide the only source of uncertainty. There are a number 
of realizations of the LISA orbits that can be used to derive expected Doppler frequencies. All 
exhibit a primary frequency of ~ 1 yr _1 with harmonics of various amplitudes. There is also 
a secular component that tends to degrade the constellation (higher Doppler shifts, larger arm 
length mismatches, etc.) as the mission progresses. In Figure [2] we plot each of the six Doppler 
parameters in (l20l) resulting from an orbital solution by Hughes (l9| that was optimized to minimize 
the average Doppler frequency in each arm. Each plot contains three traces, one for each possible 
choice of master SC. 
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B. Doppler Estimation Methods 

A number of methods have been proposed for determining the Doppler frequency. One method 
is to use the orbital ephemeris, such as the one plotted in Figure [2] to predict the Doppler. With 
periodic updates to the ephemeris from ranging data taken during normal SC down-link operations, 



the ephemeris velocities can be expected to be accurate to 3mm/s 20|. One issue is that the 
measured velocity is the projected component along the line of sight between Earth and the SC. 
Transverse velocities are not directly measured but are still constrained by the orbital model. 
Consequently one would expect that the errors in Doppler estimation could differ by a large 
amount between different inter-SC links. 

Another method for estimating the Doppler frequency is to differentiate the active ranging 
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signal that is used to determine the absolute link lengths for the TDI algorithm. Unlike ground 
tracking, this method directly measures the velocity along the inter-SC link. Ranging is expected 
to have position accuracy of ~ 1 m or better over averaging periods of ~ 1000 s |2l). This suggests 
that velocities could be measured to ~ mm/s accuracy, corresponding to Doppler frequency errors 
on the order of ~ kHz. Additional processing such as longer averaging Kalman filtering, or 



combination with an orbital model may allow for further improvements [22|. In all cases, the 
processing (including the determination of range from the pseudo-random code) would take place 
on ground. Consequently there would be some delay before the updated Doppler model could be 
uploaded to the SC. 

McKenzie, et al. jsj] proposed a simple method for determining the heterodyne frequency directly 
from the frequency meter data itself. If we consider the expression for the main science signals 
(fTTj) . all of the terms contain mean- zero stochastic processes with the exception of the Doppler 
terms. Applying a simple averaging filter to this signal can suppress the noise terms to reveal 
the heterodyne frequency. This simple algorithm relies only on information from the master SC 
and could be easily implemented on board. Table [I] gives the errors in the Doppler coefficients 
estimated by McKenzie, et al. assuming a MZ stabilized laser with a frequency noise spectrum 
given by (fT5j) and a 200 s averaging time. 



TABLE I: Doppler errors from 200 s averaging of science signal with Mach-Zehnder stabilized laser 
frequency noise given by (Tl5|) . Adapted from Table III of (9|. A * indicates the error was greater 
than the expected signal and hence the measurement is not used. 



Parameter 


^0+ 




7o+ 


7o- 


a + 


a - 


Error 


45 Hz 


0.51 Hz 


2.2Hz/s 


0.02 Hz/s 


* 


* 
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IV. ARM LOCKING SIMULATIONS 



A. Sensor Design 

As mentioned in the introduction, a number of arm locking variants have been proposed. They 
differ in the way the science signals from the two arms extending from the master SC are combined 
to form an error signal (and of course the matching controller design that completes the control 
system). In the language of (|14p . the sensor vector S differs for each design. The general goal 
in designing the sensor vector is to make the transfer function from laser frequency noise to arm 
locking sensor output as simple as possible with flat amplitude and phase responses. Ideally, 
\P(f)\ « 1 and dlP(f)/df = 0, where 

P(f) = (21) 

For example, the original proposal for single arm locking, with sensor matrix S S i ng i e — [~, 0], 
has P S ingie{f) — i sin(27r fr)e~ 2m ^ T , where r = r 12 + T%\ is the round-trip light travel time. This 
transfer function has nulls at frequencies /„ = n/r ~ 33 mHz • n, n = 1, 2, 3..., a number of which 
lie in the LISA measurement band. Since the sensor cannot measure the frequency fluctuations 
at these frequencies, the control system cannot correct for them. Furthermore, the rapid swings 
in the transfer function phase at the f n frequencies make it difficult to design a stable controller 
that extends beyond f\. More sophisticated arm locking sensors, such as the modified dual arm 
locking; sensor (MDALS)fl make a careful blend of the two science signals to generate a sensor 
with a nearly-flat transfer function in the measurement band. 

The problem of blending of multiple sensors to generate the best possible measurement of a 



state variable is a classical problem in control theory. In a previous work(23|, we applied Kalman 
filtering techniques to generate an arm locking sensor. We will briefly review this approach here. 

We begin by making a time- invariant, discrete-time linear state space model of the LISA con- 
stellation. The state vector represents the time history of the laser frequency noise over the storage 
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time in the arms, 



u i[(k-l)At] 
v 01 [{k-2)At\ 



> . 



(22) 



v i[{k-T 12 )At] 

where At is the discretization time, k is the time index and T\j = round[(rij + Tji)/At] is the 
index corresponding to the round-trip light travel time between SC± and SCj. We assume without 
loss of generality that T\i > T13. Note that the first element in the state vector represents the 
frequency delayed by a single time step as opposed to the instantaneous frequency. 

The state vector is updated using the following relations 



~£ k +i = A ~£ k + B u k + r w k , 



(23) 



A 





1 
1 
1 

••■ 



; 

1 



(24) 



B = T= [1 ■■• 0] T . (25) 

The scalar u k represents the frequency control signal applied to the laser, in this case it would 
be the output of the block labeled G\ in Figure [TJ 



Uk = VBi(kAt) Cg> G 



i- 



(26) 
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where <g) denotes convolution. Similarly, the scalar Wk represents the instantaneous frequency 
noise applied to the laser at time kAt, 



w k = v L1 (kAt). (27) 

In words, ( l23|) says that the laser control signal and the laser noise effect only the first element 
of the state vector and the remaining elements are determined through simple time delays. 

The two frequency meter outputs can be combined to form a two-element measurement vector, 
y k which is determined from the following relation, 



is A12 (kAt) 
v A13 (kAt) 



C a^fc + D Uk + H Wk + rXfc. 



C 



••• 
••■ 

D = H = 



■1 

1 
1 



;28) 



(29) 



(30) 



C is a 2 x Tyi matrix describing how the state vector couples into the frequency meter measure- 
ments. All elements of C are zero with the exception of (1, T 12 ) and (2, T 13 ), which are —1. This 
represents the fact that the frequency meter measurement includes a copy of the master laser phase 
delayed by the round-trip light travel time in the arm. The 2x1 vector represents the noise 
in each of the two frequency meter signals at the current time step. This includes the shot noise, 
clock noise, and spacecraft jitter noise contributions. The noise level can be determined from the 
noise spectra in section Hi Fl and the transfer functions to frequency meter outputs in (TTBl) . 

Equation s f[23l) and (128]) are the classical state-space description of a linear system. Kalman 



filtering 
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25] is a prescription for generating an optimal estimate of the state vector provided 



information about the system matrices and noise processes are available. In our case, the state 
matrices are determined by the arm-lengths (through the definitions of T\% and T13) and the noise 
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models that determine Wk and n The output of the Kalman filter is an estimate of the state 
vector, Xk, the first element of which corresponds to the laser frequency at time (fc — l)At. This 
element is the output of our Kalman-filter based sensor, which we refer to as an Optimal Arm 
Locking Sensor or OALS. 

In Figure [21 we show a comparison of the transfer function from laser frequency noise to sensor 
output for the single-arm, MDALS, and OALS designs. The round-trip arm lengths were chosen 
to be 33 s and 32.4 s for all three sensors for direct comparison purposes. There is no specific 
significance to the choice. The OALS sensor was computed with At = 0.1s and perfect arm- 
length knowledge was assumed for both MDALS and OALS. When compared with the single-arm 
sensor, both MDLAS and OALS exhibit much flatter responses in the LISA measurement band. 
This allows arm-locking systems based on them to achieve more uniform suppression, particularly 
near the round-trip frequencies. The OALS has less ripple than the MDALS at frequencies below 
~ 100 mHz but reaches a peak in-band ripple around 300 mHz which is similar to that of MDALS. 
Both sensors can be used to build arm locking systems that meet the LISA performance criteria. 
The OALS is optimal in the sense that it is generated using optimal control theory techniques. 
When paired with a suitable controller, we find that the net system performance is similar to that 
with achieved with the MDALS sensor, which gives us some confidence in that design. 
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FIG. 3: Transfer function from master 



aser frequency noise to sensor output for various arm locking 
sensors: Single arm sensor MDALS q|, and OALS (this work). For all cases, the round-trip light 
travel times in the two arms are 33 s and 32.4 s. 



B. Controller Design 

The second component in an arm locking system is a controller, which takes the estimate of the 
laser frequency provided by the arm locking sensor and generates a frequency tuning command 
for the laser. The design goals of the controller are to provide sufficient gain within the LISA 
measurement band to suppress the intrinsic frequency fluctuations of the master laser (TT5]) below 
the levels tolerated by TDI ([I]). As mentioned in section [TTT| care must also be taken to minimize 
the controller gain at very low (below measurement band) frequencies to mitigate laser frequency 
pulling. 

The controller is based on a classical lead-lag design. It includes a second-order lead filter at 
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the lower frequencies (break frequency at 0.05 mHz) to abate laser pulling due to uncompensated 
Doppler and Doppler derivative. It also includes a shaping filter and a single-order attenuation 
filter at 4 Hz to limit the controller action to the LISA band. Figure H] contains a Bode plot of the 
controller. 




1CT 6 KT 4 10~ 2 10° 10 : 

f[Hz] 

FIG. 4: Bode plot of arm locking controller 



C. Simulation Design 

We implemented a discrete-time simulation of arm locking as described in the preceding sections 
using the SIMULINK software package. Each arm was modeled in a manner consistent with section 
HTl The round-trip arm lengths were assumed to be t\% + T21 = 33s/c and T13 + T31 = 32.4 s/c, 
where c is the speed of light. The phase lock loops on SC2 and SC3 were assumed to be perfect 
(G*2 = G$ ^> 1) with constant frequency offsets (vm21 — 10 MHz, VM3i = 15 MHz). The Doppler 
shifts in each arm were modeled as a linearly-varying frequency with the coefficients provided 
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by the orbital model in Figure [2] at a time t = 1 yr. Doppler errors were linear in time with 
the coefficients provided in Table [I] The spectrum of intrinsic frequency fluctuations in the laser 
systems was modified from ( fl5l) to include two poles at 0.6 /iHz, limiting the total frequency 
excursion to ~ (20 MHz) over the maximum simulation period of two weeks. A two pole roll off at 
0.5 Hz was added to the spacecraft jitter noise in ([7J) to model the dynamics of the SC above the 
measurement band. 

The simulation cadence was 500 fis. System dynamics, noise generators, and the controller 
operated at this cadence. The OALS was implemented with the designed 10Hz sampling rate, 
with appropriate downsampling and upsampling filters providing the rate transitions. The OALS 
filter order was also reduced from the nominal order of 332 to 38 using balanced reduction^ |. 
This reduction provides a dramatic increase in simulation speed without changing the behavior in 
the LISA measurement band. 

V. RESULTS 

A. Component Noise Sources 

The first goal of the time-domain simulation was to verify the analytic, frequency-domain model 
of the arm locking system. Figure [5a] contains a noise decomposition of the OALS arm locking 
system derived from an analytic model. As can be seen, the overall noise in the stabilized laser is 
dominated by residual laser frequency noise, with the other noise sources being nearly four orders 
of magnitude smaller. Figure [5b] shows a similar plot obtained using the time-domain simulations. 
To obtain each curve, the simulation was run with all noise sources except the source of interest 
turned off. In all cases, the Doppler estimation errors were set to zero. The time series were then 
used to estimate a spectra. The two plots show good agreement over most of the LISA band. The 
primary differences are a broadening of the sharp spectral features near / = n/r and a roll-off 
at low frequencies in the time-domain plots. Both of these effects are consistent with spectral 
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estimation errors in the logarithmic power spectral density algorithm 27[ used to compute the 
spectra from the time series outputs. 




B. Laser Frequency Pulling 

The second goal of the time-domain simulations was to explore phenomena that are not easily 
treated analytically in the frequency domain. The laser frequency pulling described in section II III 
is an important example of such a phenomenon. We ran a simulation spanning two full weeks (the 
expected time between SC maintenance periods) with Doppler estimation errors consistent with 
our models of the errors in the averaging method. Figure [6] shows the results of this simulation. In 
the top panel, there are two curves plotted: the frequency change of the arm locked system and the 
frequency change of the intrinsic MZ noise. The first thing to notice is that the arm locked system 
drifts over approximately 20 MHz over the two week simulation period, well within the expected 
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linear tuning range of the LISA lasers. The second thing to notice is that the frequency drift in the 
arm locked system is approximately equal to the drift in the intrinsic noise. This is due to the fact 
that the arm locking loop has no effect below the LISA measurement band. The lower panel plots 
the difference of the arm locked and intrinsic frequency drifts, which gives an indication as to the 
level of additional drift generated by the arm locking system. After an initial transient decays over 
the first few days, the remaining fluctuations are less than 1 MHz. This demonstrates that this 
arm locking design does not produce any significant pulling of the master laser frequency beyond 
what is already present in the MZ stabilization system. 
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FIG. 6: Top panel: comparison of arm locked laser frequency drift with intrinsic laser frequency 
drift. Bottom panel: Difference of the two curves in the upper panel giving a rough estimate of 
the contribution to laser frequency drift from the arm locking system. 
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C. Robustness to Arm Length Errors 



Like the MDALS sensor, the OALS sensor requires some a priori knowledge of the LISA arm 
lengths. Many of the same techniques described in section IIIIBI that can potentially be used for 
estimating the Doppler frequencies can also be applied to estimate arm lengths. Which technique 
is most applicable will depend on how sensitive the performance of the arm locking system is to 
errors in the estimated arm lengths used to compute the sensor. For example, if the maximum 
tolerable error is ~ 1 m then active ranging is likely the best candidate. If, on the other hand, 
errors of ~ 10 km are tolerable, it may be possible to compute them from orbital ephemerides on 
the ground and upload new coefficients to the OALS periodically. 

To check the robustness of the OALS to errors in the arm length, we first define the mean and 
differential arm lengths assuming SC\ is the master SC, 

T m = - [Ti2 + T 21 + 7i3 + T31] , (31) 
5t = [T12 + 721 - 713 - T31] , (32) 

We then design an OALS for a specific set of nominal arm lengths, t4 0) = 32.85 s and 6t<® = 0.3 s, 
corresponding to the constellation geometry at t ~ 1.25 yrs in the orbital solution used in Figure 
|2j This sensor is used to stabilize an array of arm locking systems with different true arm lengths, 
which corresponds not only to an error in the determination of the true arm lengths, but also 
represents the situation where the spacecraft constellation has evolved in time away from the 
design point. To quantify the effect of arm length errors, we define the figure of merit 



Vpre-TDl(f) 

max- 



/ voi(f,T m ,Sr)_ 



(33) 



y (T m ,5r) = 20 log 10 

where Voi(f, T m , 8t) is the residual noise in the master laser and v pre -TDi{f) is the residual 
noise requirement specified in ([1]). \l/o measures the minimum margin in the LISA measurement 
band between the arm locking system system and the pre-TDI requirement. Where this minimum 
is positive, the frequency stabilization is guaranteed to meet performance requirements. 
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Figure [7] shows a contour plot of \l/o for this example plotted on the (r m , St) plane. The design 
point (rm,Sr^) is indicated by a white diamond. The evolution of r m and St due to the LISA 
orbit near the design is indicated by the dashed line with the grey dots indicating time intervals 
of 15 days. Figure [7] shows that positive margin exists for ~ 20 days prior to the design point 
and ~ 100 days afterwards. This is much larger than the expected intervals between maintenance 
activities, indicating that updating the OALS sensor coefficients will not drive the maintenance 
schedule of the mission. The existence of a large operating window gives us some confidence that 
the OALS design is robust enough for the LISA application. 

When arm length errors eventually do cause the system performance to degrade to the point 
where the margin is inadequate, it will be necessary to change the OALS coefficients and possibly 
which SC is designated as the master. Changing the coefficients can likely be done smoothly 
without losing lock or degrading system performance. Changing the master SC will require re 
configuring of the phase lock loops aboard all SC and will result in some down time. This should 
only be required when St for a certain arm combination becomes sufficiently small, likely 1 — 3 
times per year. 
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FIG. 7: Robustness of an example OALS-based arm locking system to errors in arm length versus 
the mean (r m ) and differential (St) arm lengths. The contours show ty Q , the minimum margin 
in any given frequency bin within the LISA measurement band. The white diamond marks the 
design point for the sensor and the dashed line shows the evolution of the arm lengths due to 
LISA orbital motion near the design point with the grey dots spaced in time by 15 days. The 
dark red color shows a positive margin between and 5 dB indicating that the performance of the 
system meets or exceeds the requirement. Evolution of the system is from left to right, indicating 
that the system has adequate performance from approximately 20 days before the design point to 
approximately 100 days after. 
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VI. CONCLUSIONS 



Arm Locking is a candidate laser frequency stabilization technique for LISA. We have used a 
time domain simulation to study a design for an arm locking system based on a Kalman filter 
optimal blended sensor and a controller which meets the frequency stability requirements for LISA 
assuming the master laser is pre-stabilized to a level of 800Hz/\/Hz. Time domain simulations 
allowed us to study transient phenomena and the performance of the stabilization system as the 
conditions of the LISA constellation evolve during the normal orbital motion, including Doppler 
shifts and imperfectly known arm lengths. The simulations indicate that it is possible to implement 
arm locking without excessive pulling of the the master laser frequency, and that the arm locking 
sensor performance is robust against errors in the absolute arm length estimates. This robustness 
allows the sensor to be periodically updated with pre-computed filter coefficients at intervals that 
are operationally reasonable. Although our time-domain simulation necessarily included specific 
sensor and controller designs, the same simulation infrastructure could be applied to study other 
candidate sensor and controller designs. 
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Appendix A: Key to Notation 



TABLE II: Partial key to notation and comparison with |!) 



Symbol 


Description 


Correspondence in [9j 


Gi(s) 


Controller transfer function on SCi 


Gi{s) 


S 


Arm locking sensor vector 


S 


— V 

v t 


total velocity of SCi 


N/A 




orbital component of SCi velocity 


N/A 


—> 
SVi 


jitter component of SCi velocity 


(6V i -6V j )-fj ij = -§ i AX ij 


Uf) 


fractional frequency fluctuations of SCi clock 


Uf) 


fjij 


unit vector from SCi to SCj 


N/A 


A 


wavelength of lasers 


A 




output of frequency meter ij 




VB1 


output of arm locking sensor 




VCij 


clock noise generated by frequency meter ij 




VDij 


Orbital Doppler shift measured by frequency meter ij 


N/A 


VElj 


Error in heterodyne model 


VDElj 


VHij 


heterodyne signal at photoreceiver ij 


N/A 


Vjij 


Spacecraft jitter Doppler shift at measurement ij 




VLi 


intrinsic frequency noise of Laser i 


§i<t>Li 


VMij 


Heterodyne model signal on frequency meter ij 


\i i = 2,3 


v M+{-) 


model of common (differential) component of heterodyne signals on SC\ 


N/A 


vo+(-) 


constant part of vm+{~) 


u o+(~) 


7o+(~) 


linear part of Vm+(—) 


7o+(-) 


a o+(-) 


quadratic part of z^v/+(-) 


«o+(-) 


VOi 


frequency output o|(Laser i 




VSij 


shot noise at photoreceiver ij 


§t4>sij 


Tij 


light travel time from SCi to SCj 


Tij 
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